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SUMAIARY 


A  surface-wave  instability  theory  applicable  to  spin¬ 
ning  discs  and  cups  is  presented  and  a  new  explanation  of  atom¬ 
ization  from  these  devices  is  proposed.  The  theory  predicts 
an  unstable  disturbance  wavelength  that  correlates  with  experi¬ 
mental  results  of  Hinze  and  Milborn  for  ligament  formation 
from  a  rotating  cup.  The  linear  analysis  is  restricted  to  flows 
where  perturbation  amplitude  is  much  less  than  film  thickness. 
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ON  ATOMIZATION  AND  LINEARIZED  FREE-SURFACE 
INSTABILITY  ON  ROTATING  BODIES 


1.0  INTRODUCTION 

Flow  with  a  free-surface  on  rotating  bodies  is  of  interest  in  the  study  of  rotary 
atomizers.  Such  devices  find  application  in  agriculture  for  dispensing  insecticides  and 
herbicides,  in  combustion  research  and  furnaces  for  spraying  fuel  in  a  finely  divided  state 
into  combustion  chambers,  and  in  spray-drying  processes  of  many  kinds. 

The  purpose  of  this  paper  is  to  develop  a  surface-wave  instability  theory  to 
explain  the  process  of  atomization  from  spinning  discs  and  cups.  The  motivation  to  do 
this  rests  on  the  proven  surface  instability  of  fluid  ligaments  and  sheets  and  in  the  hope  of 
showing  a  unified  mechanism  for  the  explanation  of  fluid  atomization.  No  attempt  is  made 
to  discuss  drop-size  distributions  resulting  from  the  atomization  process  and  no  new 
experimental  data  are  presented.  However,  attempts  at  using  existing  experimental  data 
to  verify  the  theory  are  made. 

The  steady  solution  of  Drummond'*  for  the  viscous  flow  with  a  free-surface  on 
rotating  bodies  is  used  as  a  starting  point  and  the  work  of  Lamb^’  is  used  as  the  reference 
for  the  free-surface  perturbation  analysis. 

Pedley^’  has  studied  the  stability  of  routing  flows  with  a  free  surface  where  the 
free  surface  is  cylindrical  in  form  and  can  be  either  an  "inner"  or  "outer"  boundary. 
His  linear  analysis  is  limited  to  inviscid  fluids  with  "large"  depth  in  conjunction  with 
"inertia-less"  air.  However,  he  derives  sufficient  and  sometimes  necessary  conditions 
for  stability  by  a  method  quite  different  from  that  of  Lamb.  Pedley  provides  an  extensive 
list  of  references  to  tree-surface  instability.  His  results  will  be  extended  in  this  work  to 
cover  approximately  viscous  flows  by  Lamb’s  method  and  the  two  results  will  be  shown  to 
be  compatible. 

Consider  the  flow  of  a  viscous  fluid  on  a  rotating  disc  (Fig.  1).  The  depth  of  the 
film  varies  with  radius  on  the  flat  part  of  the  disc  and  also  with  or  on  the  circular  arc 
edge.  However,  as  long  as  a  is  a  bit  less  than  90®,  then  the  film  is  nearly  of  constant 
thickness  on  the  edge.  Figure  2  shows  a  typical  thickness  variation  with  r  and  or  from 
Reference  1.  Also,  if  the  radius  of  curvature  r^  of  the  circular  arc  edge  is  small,  then 
the  radius  r  to  any  point  on  the  arc  is  nearly  constant  and  is  not  a  function  of  a.  The 
velocity  component  v,  in  the  r  direction  gets  very  small  on  the  edge,  and  the  tangential 
velocity  is  the  dominant  component.  Figures  shows  typical  values  from  Reference  1. 
The  normal  component  is  very  much  less  than  v, .  We  wish  to  investigate  surface  waves 
on  the  circular  arc  edge  under  the  assumptions  of  constant  tangential  velocity  and  con¬ 
stant  film  thickness.  At  first,  the  liquid  film  will  be  considered  inviscid  and  of  large 
depth  and  later  the  effects  of  viscosity  will  be  included. 


2.0  INATSCID  SOLUTION 

With  reference  to  Figure  4,  we  let  y  =  0  correspond  to  the  undisturbed  free 
surface  and  perturb  the  surface  an  infinitesimal  amount  r].  We  constrain  -h  <  y  <  h' 
where  h  is  the  film  thickness.  Later,  we  will  let  h’  —  w,  since  the  air  is  assumed  to  be 
unbounded  around  the  rotating  disc. 

Let  the  liquid  have  tangential  velocity  Vj  and  density  p,  and  let  the  air  have 
corresponding  parameters  V*^  and  P’.  The  surface  tension  of  the  liquid  is  T.  The  per¬ 
turbation  potentials  (f)^  of  the  liquid  and  0’,  of  the  air  must  satisfy  Laplace's  equation, 


and  the  pressure  must  satisfy  a  surface  condition  to  be  discussed  later  at  y  =  0.  In 
addition,  the  velocities  normal  to  the  free  surface  of  liquid,  air  and  the  surface  itself 
must  be  compatible.  No  flow  is  allowed  when  y  =  -h  or  h’  since  the  disc  is  solid  (y=-h) 
and  the  disturbance  potential  of  the  air  (<p\)  must  be  zero  at  the  boundary  of  the  air  (y=h’) 

1  Disturbance  Velocity  Potential 

To  satisfy  the  above  requirements,  we  assume  the  following  forms  for  the  dis¬ 
turbance  potentials  0,  and  <p\  (Lamb^'): 

(p^  =  C  cosh  [k(y+h)]  (1) 

<p\  =  C  cosh  [k(y-h»)]  (2) 

where  k  =  2n/\  and  X  is  the  wavelength  of  the  disturbance.  These  forms  satisfy  Laplace 
equation  (V^0,  =  0,  =  0),  and  the  velocity  conditions  (9</),/9y  =  0  at  y  =  -h, 

9<^»,/9y  =  0  at  y  =  h’). 

2  2  Pressure  Condition  at  the  Free  Surface 


Surface  tension  forces  allow  a  curved  surface  to  support  a  pressure  drop 
according  to  the  law 


where  R,  and  Rj  are  the  free  surface  radii  of  curvature  in  the  r,0  plane  and  axial  di¬ 
rections  respectively  (Lamb**  p.  471).  Now,  1/R,  =  1/r’  -  (d*r*/d0*)/r*  where  r’  is 
the  perturbed  free  surface  shape.  Here^  r’  =  r+q  and  rj  is  the  surface  perturbation 
relative  to  y  =  0.  From  Lamb, 


p  _  pt 


1  9^(r-h?)  _  9*(r-h?)\ 

(r-H?)*  90^  9z^  / 


(3) 


In  Equation  11  we  will  neglect  variation  of  rj  with  z;  i.  e. ,  in  the  direction  of 
the  spin  axis  (Rj  =  «)•  Expanding  Equation  3  and  retaining  only  first  order  quantities  of 
T]  and  its  derivatives  results  in 


We  assume  that  t/«9*tj/90^.  This  term  could  be  retained  but  little  accuracy  is 
gained  when  the  angular  velocity  is  large.  The  perturbation  terms  p  and  p’  are  hence 
related  by  the  perturbation  equation 

P  .  p.  +  X  ^  =  0  (4) 

r*  96^ 

2.3  Surface  Conditions  for  Velocity 


The  conditions  that  the  normal  velocity  of  the  fluids  be  the  same  as  that  of  the 
free  surface  for  small  deflections  and  that  the  tangential  velocity  is  the  only  significant 
velocity  component  leads  to; 

d<P, 


V. 
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2.4  The  Equations  of  Motion 


The  perturbation  quantities  must  satisfy  linearized  forms  of  the  equation  of 
motion  (Lamb''  ) 


P,  +  p 
P 


+  P’ 


90,  r/  1 

TT  -  [v^T  -T  -m-)  ^  \w)  j  ■  (-rn'sino!+gcosa;)T, 

dc}>\  r/  90»,y  /90’An 

—  -  h  -  -  -p-j  ^  j  J  ■  sina+gcosQ')T] 


(7) 

(8) 


where  p  and  p’  are  perturbation  pressures  relative  to  the  steady  pressures  P,  and  P’,. 
The  centrifugal  force  modifies  the  usual  force  potential  in  the  last  term  of  Equations  7 
and  8.  When  Q  is  large,  it  is  reasonable  to  neglect  g  with  respect  to  rfl^.  In  any  event, 
the  angle  a  is  nearly  90®.  Expanding  Equations  7  and  8,  neglecting  steady  terms  and 
linearizing  we  get  the  following  perturbation  equations: 


p/p 

90,  V.  8<#), 

(9) 

p»/p» 

90’,  V\  90’, 

"  ~  TT 

(10) 

Under  the  assumptions,  sina  is  taken  to  be  one. 

2.5  Surface  Shape 

We  assume  that  the  perturbation  to  the  surface  elevation  can  be  written 

Tj  =  (11) 

Note  that  9^rj/90^  is  (kr)^7j  justifying  the  approximation  in  Equation  4  for  large  k. 

2.6  Development  of  Stability  Criteria 

We  substitute  Equations  1  and  11  in  Equation  9  and  obtain 


p/p  =  [C  cosh  |k(y+h)|  ji(a-kV^)|  +  rn*a]e'*‘’'''‘'®’  (12) 


Similarly,  Equations  2  and  11  substituted  into  Equation  10  yields 


p’/p’  =  [C»  cosh  |k(y-h»)}  {i((T-kV»^)}  +  rn2a]e‘''''-'''«' 


(13) 
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We  apply  the  pressure  condition,  substitute  Equations  12  aiid  13  into  Equation  4,  set 
y  =  0.  retain  only  perturbations  and  simplify.  The  result  is: 


pc  cosh  kh(i)((r-kVj)  -  P'C»  cosh  kh'(i)(cr-kV',)  =  a[Tk^  -  (P-P')  rfi^] 
The  sui'face  conditions  (Eqs.  5  and  6)  yield: 

C  =  “  ni(a-kVj)/(l<  sinh  kh) 

C»  =  +  ai(cr-kV»,)/(k  sinh  kh') 

Substituting  Equations  15  and  16  in  Equation  14  results  in: 
p  coth  kh^^  -  V,y  +  P'  coth  kh' =  Tk  -  (P-P')  rfi^/k 
We  solve  for  (cr/'k)  in  Equation  17  and  let  S  =  P'  coth  kh'/p  coih  kh  and  get 


(14) 

(15) 

(16) 

(17) 


V,  +  SV'^ 

1  +  s 


1  +  £ 


(18) 


-(Vt-W,)'  -  (1+S)(P/P')  [i|-  (I-P'/P)  - 


The  wave  speed  is  or/k. 

Examination  of  Equation  18  shows  that  a  will  be  complex  (y±i6)  if  the  quantity 
f(k)  under  the  root  sign  is  negative.  From  Equation  18, 


f(k)  =  -  (V,-W,)'  -  (1+S)(P/P')  (I-P'/P)  - 

If  f(k)  is  negative,  then  the  perturbation  amplitvide  ri  becomes 


(19) 


rj  =  ae'"''*''’’'''  e 


61 


which  is  divergent.  Hence,  if  f(k)  from  Equation  19  is  less  than  zero,  the  wave  is 
unstable.  Setting  f(k)  equal  to  zero  defines  a  boundary  for  stability.  In  the  present 
application,  the  relative  velocity  between  the  liquid  and  tlie  stationary  air  (V, -W,)  is 
just  rS2.  Further,  h'  <»  and  the  factor  S  reduces  to  (p'/p)  tanh  kh  which  is  «1.  The 
density  ratio  p'/p  is  also  very  small  for  air-liquid  combinations.  Hence,  we  approximate 
Equation  19  by 

f(k)  =.  -  (rn)'  -  (p/p')  0  (20) 

and  obtain  an  approximate  stability  criterion.  Equation  20  is  quadratic  in  k.  We  define 


A*  = 


27rr 


(kr) 


(21) 


as  a  nondimensional  disturbance  wavelength.  Substitution  of  Equation  21  in  Equation  20 
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leads  to  the  following  expression  for  f(A*) : 

f(X*)  +  (P»/P)X*  -  1/W  >  0  (22) 


where  W  is  Weber's  Number  defined  by 

W  =  (23) 

There  is  one  real  positive  root  to  Equation  22  denoted  by  A*, ,  Figure  5  illus¬ 
trates  the  range  of  A*  for  stable  or  unstable  waves. 

Solving  Equation  22  for  A*,  leads  to  the  following  limit  on  A*  for  unstable 
surface  waves: 


A*  >  A* 


P'/P 


/P)  wj 


(24) 


If  A*  <  A*  an  unstable  surface  wave  results. 

Pedley^’  neglected  the  air  inertia  (P*  -  0)  in  his  analysis  but  he  did  not  make 
the  assumption  17  «  (Eq.  4).  He  quotes  the  result  of  Hocking^’  for  the  necessary  and 
sufficient  conditions  for  stability  as 

5^  (5^  ^  1) 

<  — when  Tj  « 


A* 


90 


(25) 


where  the  notation  has  been  made  to  agree  with  this  work.  Setting  p'  =  0  in  Equation  22 
leads  to  A*^  >  1/W  for  instability,  which  is  the  same  as  Hocking's  result.  Since  the 
maximum  value  of  A*  is  one,  and  since  there  must  be  an  integral  number  of  wavelengths 

in  order  for  a  periodic  surface  wave  to  exist,  it  is  easily  seen  that  JW  =  1,2,3, . n 

for  this  special  case.  More  generally,  we  can  include  the  air  inertia  (p'  4  0)  and  obtain 
from  Equation  22 


W  > 


n‘ 


1  +  (PVP)n 


n  =  1,2,. 


(26) 


3  ^  77 

for  instability.  Properly,  n  >  10  for  the  approximation  rj «  — Now  p'/p  is  of 

j  90 

order  10'^ for  air-liquid  interfaces,  so  that  for  n  <  10  ,  W  2:  n’  is  a  good  approximation 

to  the  instability  criterion  while  for  n  5:  lo\  W  2:  ^  is  better.  For  W  sufficiently 

high,  a  wide  range  of  disturbance  wavelengths  lead  to  an  unstable  free  surface.  Squire^’ 
has  shown  that  a  sheet  of  liquid  moving  in  still  air  is  unstable  provided  the  Weber  Number 
(based  on  one  half  the  sheet  thickness)  Is  greater  than  one.  He  quotes  fair  agreement 
with  cApei  iiViCnt. 
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2.  7  Effect  of  Parameters  on  Surface  Stability 

The  surface  tension  enters  the  problem  only  in  the  Weber  Number  W.  As 
T  —  0,  W  —  oe  and  the  stability  criterion  of  Equation  22  reduces  to 

f(X*)  =  +  P'/P)  >  0  (27) 

This  Is  satisfied  by  all  positive  A*  and  shows  the  stabilizinjr  effect  of  surface  tension  in 
the  sense  that  a  smaller  range  of  values  of  A*  result  when  T  0. 

The  result  of  Equation  27  is  reproduced  by  large  values  of  r  ci  n.  For  a 
given  disc,  the  higher  the  angular  velocity  the  wider  the  range  of  A*  for  unstable  waves. 
Hence,  fl  is  a  "destabilizing"  variable. 

The  density  ratio  p’/p  can  vary  by  choice  of  liquid  (p)  and  its  temperature  as 
well  as  by  varying  ambient  conditions  of  air  temperature  and  pressure  (p') .  The  addition 
of  P’/p  acts  in  a  destabilizing  sense  Increasing  p  acts  in  a  stabilizing  sense  while  the 
converse  is  true  for  p’. 

2.8  Maximum  Instability  with  "Large"  Film  Thickness 

Some  disturbance  wavelength  in  the  unstable  region  maximizes  the  imaginary 
part  of  a,  leading  to  a  "fastest  growing"  disturbance.  From  Equations  18  and  21  we  obtain 
the  following  expression  for  at 


^  Fr  *  ypvWanh  kh  /g(A^ 

where 

g(A*)  =  +  v\+  iP/p')v\\*  -  (P/P')  ^ 


(28) 

(29) 


The  hyperbolic  tangent  is  nearly  1  when  kh  is  greater  than  about  1.5.  In  this  sense,  we 
let  kh  be  "large",  and  let  tanh  kh  =  1  in  Equation  28.  The  result  for  the  im.'g'nary  part 
of  cr  is: 


Im((T)  -  i 

V  X* 

The  condition  for  maximum  instability  is; 

d[lm(tT)]  _  „ 
dA+  '  ” 


The  operation  indicated  by  Equation  32  is  applied  to  Equation  29  which  results  in; 


(30) 


(31) 


(32) 


A*^  +  2(pVp)A*  -  r.  0 


(33) 


for  finite  A*. 
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Equation  33  is  quadratic  in  A,*  and  it  has  one  root  (^’cru)  within  the  unstable  region: 


X* 


CRII 


(pVp) 


1  + 


3 

(P’/P)^  W 


(34) 


Disturbances  with  wavelength  equal  to  are  those  which  grow  the  fastest.  Substi¬ 

tuting  Equation  34  in  Equation  28  leads  to  the  corresponding  complex  value  of  a.  Note 
that  1A*crit  must  be  an  integer  leading  to  the  conclusion  that  very  often  there  is  no 
physically  possible  solution  to  Equation  34.  Hence  any  X*  in  the  unstable  range  is  equally 
likely.  The  higher  the  Weber  Number,  the  more  likely  a  solution  to  Equation  34  becomes. 


Integral  values  of  (X*)"'  and  (X*crit)*’  were  chesen  and  the  corresponding  Weber 
Number  for  P  =  1.65  slugs/ft^  and  T/p  =  .0015  fi  Vsee*  was  computed  from  Equations  24 
and  33.  The  results  are  shown  in  Figure  6.  Properly,  the  points  should  not  be  joined  by 
a  smooth  curve. 


For  many  liquids,  T/p  does  not  vary  outside  the  range  0.  5  x  10'^  to  2.5  x  10'\ 
Below  a  certain  speed,  the  waves  are  stable.  For  example,  for  (X*)"’  =  10,  the  waves 
are  stable  for  r^n^  <  .148,  (T/p  =  .0015).  The  smaller  the  radius,  the  higher  this 
minimum,  speed  becomes.  It  is  seen  that  X*crit  is-  close  to  X*,.  For  W  <  10^,  log  (X*)"' 
is  approximately  a  linear  function  of  log  W  with  a  slope  of  .52.  Hence,  for  inviscid 
fluids,  (X=|)''  =  W  +  constant. 


3.0  VISCOUS  SOLUTION  FOR  ’’LARGE**  DEPTH 

Lamb^’  (p.  627)  considers  surface  waves  of  a  stationary  viscous  fluid.  He 
concludes  that  the  wave  velocity  is  very  nearly  the  same  as  for  an  inviscid  solution  but 
that  the  amplitude  of  the  wave  is  attenuated  approximately  by  a  factor  e'2»'k*t.  The  par¬ 
ameter  V  is  the  kinematic  viscosity  of  the  liquid.  Harrison'^*  has  analyzed  the  viscous 
problem  more  accurately  and  we  interpret  Lamb’s  result  as  applying  when  p*p’«pt'.  The 
restriction  is  easily  met  for  air-liquid  interfaces  even  if  v^>v.  Further,  the  flow  irro- 
tationality  assumption  becomes  more  approximate  the  nearer  one  gets  to  the  disc  surface 
since  viscous  forces  increase  as  y  -*  -h.  Therefore,  we  add  a  further  condition  that  the 
results  will  be  accurate  only  when  h>2X.  In  this  event,  the  factor  S  from  Equation  18 
is  just  P’/p,  and  we  restrict  ourselves  to  large  depth  solutions. 

With  the  above  comments  in  mind,  we  replace  Equation  11  by 


n 


ae-2>'i>^*e 


(35) 


and  retain  the  wave  velocity  a/k  from  Equation  18.  A  slight  revision  of  the  Instability 
criterion  is  necessary.  As  before,  let  a  =  y  :t  i6  and  substitute  it  into  Equation  35  and 
obtain 


(36) 


From  Equation  36,  77  will  be  divergent  if 


m(k)  =  -2vk^  +  6  >  0  (37) 

Using  Equation  28  to  evaluate  6,  we  have 

-2vk^  +  k  JpVp  Jg(k)  >  0  (38) 

for  an  instability  criterion  where  g(k)  is  given  by  Equation  29.  As  before,  we  substitute 
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A*  =  (kr)"'  and  solve  for  A*: 


h(A*)  =  A*^  +  (p*/p)A**  -  (1/W)A*  -  4/Re^  >  0  (39) 


One  positive  real  root  of  Equation  39  always  exists  and,  defines  a  stability  boundary.  Two 
negative  real  roots  always  exist.  Figure  7  shows  the  stability  boundary  as  a  function  of 
llevnolds  Number  for  various  values  of  Weber  Number.  Values  of  p*  and  p  are  the  same 
as  those  for  Figure  6. 

All  the  curves  coalesce  at  low  Reynold’s  Number  which  shows  that  the  stability 
boundary  for  very  viscous  fluids  is  insensitive  to  surface  tension. 

Each  stability  boundary  in  Figure  7  tends  to  rise  linearly  and  then  roll  off  to 
the  inviscid  assymptote  as  Reynold’s  Number  increases.  The  roll-off  point  is  a  function 
of  Weber  Number.  The  linear  portion  can  be  approximated  by 

(A*)"'  =  Re  ^  +  C 

When  viscosity  is  zero  (Re  —  <»)  then  Equation  39  yields  the  same  roots  as  its 
inviscid  counterpart  Equation  23.  The  addition  of  viscosity  results  in  a  smaller  region 
of  unstable  wavelengths  and  hence  is  said  to  exert  a  "stabilizing"  influence  on  the  free 
surface. 


The  wave  velocity  of  the  unstable  viscous  fluid  is  identical  with  the  inviscid  and 
equals  the  tangential  velocity  .  The  wave  is  stationary  with  respect  to  the  disc. 


3. 1  Effect  of  Flow  Rate 


For  a  given  fluid,  disc  and  speed,  the  value  of  h  depends  only  on  the  flow  rate. 
An  approximate  equation  for  h”  is: 


h 


i  [E  /  3Q  \ 

Ui-r’  JnFr/ 


After  some  manipulation,  the  most  useful  form  for  h  is: 


h  = 


1/3 


where  a*  is  the  circumference  of  the  rotating  body  and  r  is  the  local  slope.  It  is  seen 
that  h  varies  as  Q’'"*. 


As  O'  7r/2  r  —  0  leading  to  a  large  film  thickness  but  h  is  nearly  constant 
for  a  range  which  ends  quite  close  to  a  =  n-/2.  It  is  easily  seen  that  a  very  small  film 
thickness  can  result  from  high  values  of  V,  and  Re  or  low  values  of  Q.  Hence,  small 
values  of  A*  are  physically  possible  that  are  well  within  the  limitations  that  A<2h. 


3.2  Atomization  and  Instability 

The  above  theory  on  surface  instability  is  applied  in  the  formulation  of  a  new 
theory  on  the  mechanism  of  atomization  from  rotating  bodies.  Consider  the  sketch  in 
Figure  8.  When  A*>A*,  the  free  surface  is  unstable,  rj  tends  to  h  as  time  progresses 
and  a  volume  of  liquid  equal  to  A’ h  is  trapped  in  the  shaded  area.  Since  rj  is  periodic  in 
0  (period  =  A*),  the  number  of  "sites"  of  formation  is  constant  around  the  body,  and  the 
amount  of  flow  to  any  site  is  A*Q  ft  Vsec.  Macfarlane  and  Colboum'”  have  verified  this 
exper  i  menta  lly . 
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Using  Equations  18  and  35,  we  can  write 


eMx-a 

(40) 

(41) 

If  n  is  large,  then  the  perturbation  grows  rapidly.  The  phase  shift  between  rj 
and  T)  is  about  90°  because  S2  and  n  are  large  while  \*  is  small.  The  result  is  a  large 
outward  velocity  over  the  rising  portion  of  the  wave  and  a  corresponding  large  inward 
velocity  over  the  falling  portion  of  the  wave.  Hence,  the  trapped  volume  will  leave  the 
rotating  body  normal  to  the  free  surface  with  a  considerable  amount  of  internal  rotation. 
Three  tvpes  of  atomization  have  been  experimentally  observed*  from  rotating  discs  and 
cups: 


1)  direct  drop  formation  at  low  Q, 

2)  ligament  formation  at  moderate  Q, 

3)  sheet  formation  at  high  Q. 

Ligaments  decay  into  drops  by  Rayleigh-type  surface  instability  of  rotationally 
symmetric  columns  of  fluids,  while  sheets  decay  into  ligaments  by  a  surface  instability 
(Hagerty  and  Shea®’  or  Dombrowski  and  Johns’’)  and  subsequently  into  drops  by  Rayleigh 
instability.  The  flow  rate  Q  determines  which  type  of  atomization  results.  Now,  if  A^h, 
the  volume  of  trapped  fluid,  is  greater  than  \*QAt,  the  flow  to  the  site,  then  type  1  atom¬ 
ization  results  in  an  intermittent  manner  because  it  takes  a  finite  time  to  replenish  the 
fluid  and  rebuild  the  fluid  layer  for  a  subsequent  unstable  ejection  of  fluid.  However,  if 
A^h  is  about  equal  to  A*QAt,  then  continuous  ejection  of  ligaments  must  occur  in  order  to 
satisfy  the  continuous  flow.  Obviously,  the  third  condition  of  sheet  ejection  results  when 
A*QAt  is  greater  than  A^h.  Probably  At  will  have  to  be  studied  experimentally,  but  it 
probably  is  related  to  n(A*)  (Eq.  40). 

If  the  instability  criterion  is  not  met,  then  the  linear  approximation  predicts  a 
stable  flow  which  separates  as  a  sheet  from  the  rotating  body  independent  of  the  flow  rate. 
If  the  flow  rate  is  low,  then  the  film  is  very  thin  resulting  in  a  very  short  sheet  which  very 
quickly  decays  to  drops. 

Historically,  the  accepted  mechanism  of  drop  formation  from  the  edge  of  a  disc 
has  been  to  equate  the  mass  of  the  drop  times  its  centripetal  acceleration  with  the  surface 
tension  force  based  on  the  drop  circumference.  The  correlation  function  due  to  Walton 
and  Prewett'®’  is 


dn 


constant 


dQ 


constant 


(42) 


(42) 


where  d  is  the  drop  diameter.  The  constant  in  Eq  -ation  42  did  not  correlate  well  using 
experimental  values  of  parameters  on  the  lefthand  side. 


The  only  steady  surface  tension  effects  are  the  rise  of  fluid  in  a  capillary  tube 
and  the  ability  of  the  curved  free  surface  to  support  a  steady  pressure  difference.  The 
former  requires  fluid  contact  on  both  sides  of  a  tube  normal  to  the  free  surface  while  the 
latter  is  subtracted  out  exactly  from  the  perturbation  equations.  Therefore  surface 
tension  forces  arise  only  from  altering  the  curvature  of  the  steady  free  surface  due  to 
frce-surfacc  perturbations.  Only  when  the  surface  Is  "slretchod"  should  surface  forces 
be  considered.  Hence  it  is  unsatisfying  to  equate  surface  tension  forces  whose  origin  is 
in  unsteady  phenomena  with  centripetal  accelerations  which  are  steady  when  n  is  constant 
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Rayleigh  instability  of  fluid  ligaments  arises  from  a  surface  instability  and  now  a  consist¬ 
ent  drop  formation  process  from  rotating  discs  based  on  surface  instability  is  proposed. 
Hinze  and  Milbom^’  applied  an  instability  approach  in  attempting  to  explain  and  correlate 
their  experimental  results.  However,  their  point  of  view  was  qualitative  and  was  used 
only  to  collect  variables  in  functional  groups  for  experimental  determination  of  a  constant 
and  an  exponent.  Pedley^’  hints  at  surface  instability  being  involved  in  atomization  of 
swirling  liquid  jets. 

As  the  amplitude  of  the  surface  wave  increases,  the  linearized  form  of  the 
equations  presented  here  becomes  less  accurate  and  the  of  the  acihjal  wave  may  differ 
from  the  theoretically  predicted  value.  However,  if  a  non-linear  system  is  unstable  to 
an  infinitesimal  disturbance,  so  is  its  linearized  approximation.  Hence,  the  extension  of 
the  proposed  drop  formation  process  to  real  situations  is  immediate. 

3.3  Comparison  of  Theory  with  Available  Experimental  Data 

Experimental  work  which  supplies  sufficient  data  to  calculate  is  sparse. 
Either  the  fluid  or  the  operating  conditions  are  not  completely  documented.  However,  an 
exception  is  the  paper  by  Hinze  and  Milborn^’,  although  even  their  method  of  presentation 
leaves  doubt  as  to  the  values  of  all  the  parameters  for  any  specific  data  point.  Their 
main  contribution  to  the  present  comparison  is  a  semi-empirical  relation  for  the  number 
of  ligaments  shed  from  a  rotating  cup.  They  present  the  relation 


for  the  number  of  ligaments.  K  and  n  were  found  to  be  .215  and  4/5  from  experiment. 
A  reasonably  large  amount  of  scatter  (*25%)  existed.  In  the  notation  of  this  paper  and 
noting  z  aquals  1/A*  we  get  (from  Eq.  43); 

A*  =  l/(.577Re'^w''‘‘)  (44) 


In  arriving  at  this  equation,  Hinze  and  Milborn  allowed  Q  to  vary  provided  that 
the  flow  rate  was  not  high  enough  to  produce  a  sheet  and  yet  high  enough  to  exclude  direct 
drop  formation.  In  addition,  data  for  three  different  liquids  were  presented. 

They  noted  a  ’'liquid  torus"  formed  at  the  lip  of  the  cup.  This  is  a  reasonable 
situation  since  the  steady  radial  velocity  of  the  liquid  varies  from  zero  on  the  cup  surface 
to  a  maximum  at  the  free  surface  and  the  inward  directed  normal  velocity  has  the  same 
behaviour  resulting  in  a  tendency  for  the  liquid  to  flow  around  the  edge  of  the  cup  lip. 
Hence,  the  approximation  of  or  in  the  stability  equations  being  equal  to  90“  is  appropriate. 

The  results  of  the  stability  analysis  using  parameter  values  within  the  experi¬ 
mental  range  are  compared  with  the  experimental  results  of  the  three  liquids  as  well  as 
the  overall  correlation  function  (Eq.  44)  on  Figures  9  to  11.  The  viscosity  of  the  liquids 
was  the  major  variation  between  them. 

On  Figure  9,  the  experimental  data  for  the  least  viscous  liquid  virtually  coin¬ 
cided  with  the  viscous  stability  boundary,  as  did  the  overall  correlation  function.  The 
comparison  is  very  good.  Increasing  viscositj'  (Fig.  10)  showed  the  overall  correlation 
function  close  to  the  viscous  stability  boundary  but  the  particular  liquid  curve  showed 
some  discrepancy.  The  highest  viscosity  liquid  (Fig.  11)  was  within  the  stable  region  as 
predicted  by  the  viscous  boundary,  but  was  within  the  unstable  region  from  the  inviscid 
boundary.  Apparently  non-linear  and  rotational  flow  effects  have  become  important.  In 
general,  the  comparison  is  favourable  considering  the  comments  on  scatter  and  variable 
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Q  already  made.  This  tends  to  substantiate  the  theory  of  this  report.  When  the  theory 
predicts  instability,  Hinze  and  Milbom's  data  show  ligaments.  The  only  area  where 
discrepancies  occur  is  for  very  viscous  liquids. 

A  further  check  on  experimental  results  was  made  from  a  photograph  by 
Macfarlane  and  Colboum"’  showing  atomization  by  direct  drop  formation  from  the  edge 
of  a  spinning  disc.  The  picture  was  a  shadowgraph  taken  looking  vertically  down  on  a 
13.4°  sector  of  a  one- inch  diameter  disc.  The  regularity  of  drop  formation  sites  was 
evident  ar.d  23  were  counted  at  the  rim  of  the  sector-  The  angular  velocitj'  was  313  revo¬ 
lutions  per  second.  Regarding  each  drop  site  as  a  wavelength,  it  was  easy  to  calculate 
X*  =  3.2x  10'^.  The  fluid  data  were  not  specifitid  but  by  using  T/p  of  1.2  x  10'^,  P’/P 
of  1.4  X  10'^,  and  =  10'^  (the  fluid  was  a  light  fuel  oil)  the  inviscid  and  viscous  boundary 
calculations  for  A*  were  1.5  x  10'^  and  2.2  x  10'^  respectively.  Therefore,  the  theory 
did  predict  an  unstable  surface  wave  and  the  proposed  atomization  model  was  not  contra¬ 
dicted. 

However,  Stuart'^*,  in  his  review  article  on  non-linear  effects  in  hydrodynamic 
stability,  points  out  that  Poisseuille  Flow  can  be  unstable  due  to  non-linear  effects  when 
the  flow  is  predicted  to  be  stable  by  linearized  theory.  There  may  be  an  energy  transfer 
from  the  perturbed  motion  to  the  mean  motion  by  Reynold’s  Stresses  which  distort  the 
mean  flow.  Rosenhead’^’  states  that  non-linear  terms  do  not  act  to  stabilize  a  flow  and 
that  plane  Couette  Flow  is  unstable  for  finite  disturbances  even  though  it  is  stable  for 
infinitesimal  disturbances.  The  linear  theory  is  thus  considered  to  be  somewhat  restric¬ 
tive  since,  it  does  not  account  for  any  energy  transfers  between  the  mean  and  perturbed 
motions. 


Lamb’’  (p.  417)  discusses  waves  of  finite  amplitude.  The  free  surface  of  a 
quiescent  liquid  of  infinite  depth  takes  on  an  ’’approximately  trochoidal  form"  which 
sharpens  the  crests  and  flattens  the  troughs  of  the  sine  wave  predicted  by  linear  theory. 
This  is  exactly  what  the  pictures  of  Hinze  and  Milborn  show  for  ligament  formation,  i.  e. 
the  space  between  the  ligaments  is  greater  than  a  ligament  diameter. 


The  linearized  theory  is  valid  only  when  a/A  is  small  (Lamb”  p.  417).  In  other 
words,  the  amplitude  of  the  surface  wave  must  be  much  less  than  its  wavelength.  For 
irrotationality,  h  must  be  greater  than  A/2.  Therefore,  the  wave  amplitude  is  restricted 
to  being  much  less  than  h.  For  atomization,  the  amplitude  must  equal  the  film  thickness 
which  leads  to  the  conclusion  that  the  linearized  theory  can  only  be  used  to  test  for  incip¬ 
ient  instability. 


4.0  CONCLUSION 

A  surface-wave  instability  theory  has  been  presented  to  attempt  an  explanation 
of  drop  and  ligament  formation  from  edges  of  rotating  discs  and  cups.  Sheet  formation 
has  already  been  explained  by  the  separation  of  the  fluid  from  the  edge  of  the  rotating  disc 
by  the  steady  flow  theorj'  (Drummond”).  Pedley’s”  results  have  been  extended.  Some 
correlation  between  theory  and  experimental  results  of  Hinze  and  Milborn”  and  Macfarlane 
and  Colbourn'”  have  been  shown  which  tends  to  substantiate  the  theory.  However,  the 
linear  theory  may  be  too  restrictive  and  non-linear  as  well  as  rotational  flow  effects  may 
be  large  during  the  actual  fluid  shedding  from  the  disc  edge. 


Atomization  from  sheets  and  isolated  ligaments  has  already  been  explained  by 
free-surface  instability,  and  now  atomization  from  rotating  cups  and  discs  has  been  put  in 
the  same  framework.  Some  experimental  work  is  needed  to  provide  data  on  edge  waves 
and  their  breakdown. 
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